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I. INTRODUCTION AND MOTIVATION 

Helicity-dependent parton densities (PDFs) contain 
the information to what extent quarks and gluons with a 
given light-cone momentum fraction x tend to have their 
spins aligned with the spin direction of a nucleon in a 
helicity eigenstate. The most precise knowledge about 
these non-perturbative quantities, along with estimates 
of their uncertainties, is gathered from comprehensive 
global QCD analyses [l|, Q to all available data taken 
in spin-dependent deep- inelastic scattering (DIS), with 
and without additional identified hadrons in the final- 
state, and proton-proton collisions. Extractions of helic- 
ity PDFs are based on the assumption that they factorize 
from calculable short-distance partonic scattering cross 
sections, which is expected to be a good approximation 
for processes characterized by a sufficiently large momen- 
tum scale above, say, about 1 -j- 2 GeV. Current analyses 
[l|-[H are performed consistently at next-to-leading order 
(NLO) accuracy of QCD for both the scale evolution of 
helicity PDFs and the relevant hard scattering cross sec- 
tions. 

Apart from being essential for a comprehensive under- 
standing of the partonic structure of hadronic matter, 
helicity PDFs draw much their relevance from their re- 
lation to one of the most fundamental and basic but yet 
not satisfactorily answered questions in hadronic physics, 
namely how the spin of a nucleon is composed of the spins 
and orbital angular momenta of quarks and gluons. The 
integrals of helicity PDFs over all momentum fractions x 
(first moments) at a resolution scale fi, 
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provide information about the contribution of a given 
parton flavor / to the spin of the nucleon. There are 
well-known subtleties related to the decomposition of the 
proton spin in QCD Q, for instance, Ag(fi) has the inter- 
pretation of the gluon spin contribution only in light-cone 
gauge, which is closely tied, however, to the QCD im- 
proved parton model, but otherwise is a non-local opera- 
tor. The contribution of the quark and antiquark spins to 
the nucleon spin, as summarized by the first moment of 
the flavor singlet combination AS(/i) = S/= g s A/(/0, 
is related to a gauge-invariant operator. Although orbital 
angular momenta of quarks and gluons have to be present 
in the scale evolution of longitudinally polarized quarks 
and gluons to obey angular momentum conservation in 
collinear 1 — > 2 parton splittings @, their contributions 
cannot be quantified from the experimental probes con- 
straining helicity PDFs. To this end, information about 
correlations between the nucleon's spin and transverse 
degrees of freedom of quarks and gluons has to be ac- 
quired; see, e.g., Refs. 0, @| for details. 

The relevance of helicity PDFs, and spin physics in 
general, is also reflected in more than a dozen vigor- 
ous experimental programs in the past twenty-five years, 
matched by tremendous advancements in the under- 
standing and development of the underlying theoretical 
framework. The most recent global analyses [lH4( con- 
firm early findings that the total quark spin contribution 
AE(/x) is significantly smaller than expectations from 
naive quark models even within still sizable uncertainties 
from extrapolations to the unmeasured small x region in 
Eq. ([1} . A potentially very large gluon spin contribution 
A<7(/x) > 2 (in units of H), initially thought to be a vi- 
able way to account for the "missing" nucleon spin [8| , is 
now strongly disfavored by measurements of jet and pion 
yields at BNL-RHIC (Sj, which, on the contrary, prefer 
a rather small Aglxju) in the range 0.05 < x < 0.2 
at fi ~ 5 GeV |l|, LL0( ■ Results for charm and hadron 
production in polarized lepton-nucleon scattering [ll[ are 
consistent with Ag at x ~ 0.1. However, Ag(x,fi) 
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remains to be completely unconstrained at x < 0.01 due 
to the lack of data and, depending on which functional 
form one assumes for the extrapolation to small values 
of x in Eq. (JlJ , sizable gluon spin contributions of up to 
Ag([i) ~ 1 are still conceivable QUI]- Theoretical argu- 
ments based on the color coherence of gluon couplings 
suggest that Ag(x, /i) ~ xg(x, ^i) in the limit x — > 
at some low but otherwise unspecified bound-state like 
scale fi [llj]. Perturbative evolution to larger scales [l3[ 
will change this small x behavior though. While helic- 
ity PDFs will not exhibit the strong rise of unpolarized 
PDFs driven by the 1 jx singularity in the evolution ker- 
nel, their actual small x behavior remains unconstrained 
by present data. 

The surprisingly small, perhaps even positive 
strangeness helicity PDF, as determined from semi- 
inclusive DIS (SIDIS) data with identified charged 
kaons in the broad range 0.005 < x < 0.5 [3, |2. \4\, 
has triggered quite some discussions recently [1J, Il5 |. 
If SU(3) flavor symmetry is approximately valid, one 
expects a significantly negative first moment of about 
As(/i) + As(/z) ~ —0.1 by utilizing the experimentally 
well determined hyperon decay constants F and D and 
the value for AS(/i) extracted from fits to polarized DIS 
data. Recent SIDIS data from COMPASS pj] exhibit a 
weak trend for As(x,/i) turning negative somewhere in 
the region 0.001 < x < 0.01 [la], and acquiring a large 
negative moment As(/z) in accordance with SU(3) sym- 
metry is still possible if As(x,[/,) is large and negative 
in the currently unmeasured small x region. However, 
lattice QCD results for As (/it) + As(/i) [l7| and compu- 
tations of SU(3) breaking effects in axial current matrix 
elements [l8[ point towards sizable violations of SU(3) 
flavor symmetry, perhaps even consistent with a vanish- 
ing total strangeness polarization As(fi) + As(/i) ~ 0. 
This might be explained by significant chiral corrections 
as was estimated in [lj| within the framework of the 
"cloudy bag model". To complicate things further, all 
current extractions of strangeness helicity PDFs from 
SIDIS data exhibit a significant dependence on the 
choice of strangeness-to-kaon fragmentation functions 
(FFs) [l6j], which needs to be scrutinized further. Forth- 
coming data from i3-factories 20], DIS multiplicities 
[2lj . RHIC, and the LHC are likely to considerably 
improve our knowledge of FFs soon. Another surpri sing 
outcome of analyses of recent COMPASS data [la, LL6[ 
was the diminishing evidence for a sizable asymmetry 
in the light quark sea, i.e., Au(x,fi) — Ad(x,fi) ^ 0; 
uncertainties are still large though. Within the large- Ac- 
limit of QCD as incorporated in, e.g., the chiral quark 
soliton model [22| one expects an SU(2) breaking which 
is at least as large as what has been already observed 
for unpolarized PDFs 23] but with the sign reversed. 

Clearly, despite the impressive progress made both ex- 
perimentally and theoretically many fundamental ques- 
tions related to the proton's helicity structure, including 
a quantitative understanding of the decomposition of the 
proton's spin still remain unanswered. The discussions 



above exemplify the need for measurements which are 
sensitive to smaller values of x than accessible with past, 
present, and upcoming polarized DIS fixed-target experi- 
ments or with high transverse momentum probes at BNL- 
RHIC. An accurate determination of the first moments 
AE(/x) and Ag(fi) entering the proton's spin sum rule 
or elucidating the flavor dependence of helicity PDFs to 
quantify, e.g., a potential SU(3) symmetry breaking in 
the light quark sea cannot be achieved without consider- 
ably enlarging the kinematic coverage of spin-dependent 
data in the future. All the required measurements to 
address and answer these questions related to the small 
x regime are unique to a polarized, high energy lepton- 
nucleon collider such as the proposed electron-ion collider 
(EIC) project 0. 

In the remainder of this paper we will demonstrate 
quantitatively how polarized DIS and SIDIS measure- 
ments at an EIC will improve our knowledge of helicity 
quark and gluon densities and their contributions to the 
spin of the nucleon. Our assessment of the impact an 
EIC is expected to have on the determination of helicity 
PDFs is based on a series of global QCD analyses at NLO 
accuracy performed with realistic sets of pseudo-data for 
inclusive and semi-inclusive DIS measurements at an EIC 
at various conceivable center-of-mass system (c.m.s.) en- 
ergies. Uncertainties of PDFs are estimated with both 
the robust Lagrange multiplier method [24[ as well as 
within the Hessian approximation [25[. We will briefly 
touch upon related inclusive measurements such as the 
Bjorken sum rule, the charm contribution to the polar- 
ized structure function gi(x,fj,), and novel electroweak 
probes in spin-dependent DIS. 

In the next Section we will describe how the projected 
data for DIS and SIDIS are generated, their kinematic 
coverage, and what kind of cuts have been imposed. We 
also outline our method to quantify the impact of the EIC 
data on determinations of helicity PDFs. Section III is 
devoted to detailed discussions of the results of the global 
analyses performed with projected EIC data. Other op- 
portunities related to the nucleon's helicity structure will 
be briefly discussed in Sec. IV. The main results will be 
summarized in Section V. 



II. KINEMATICS, STRATEGY AND 
FRAMEWORK 

An EIC will most likely be realized in at least two 
stages with increasing c.m.s. energies [7(. To assess the 
impact of a future EIC in determining helicity PDFs we 
will consider two sets of energies conceivable with the 
first stage of the eRHIC option of an EIC (2(| which is 
based on colliding an E e = 5 GeV electron beam with the 
existing RHIC proton beam of E p = 100-250 GeV. Sim- 
ulations based on pseudo-data generated with an elec- 
tron energy of 20 GeV are used to estimate the impact 
of a later stage of an EIC. The resulting c.m.s. ener- 
gies y/s and corresponding lowest accessible values of 
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x = Q 2 1 (sy) for two different values of momentum trans- 
fer Q 2 are summarized in Tab. HI assuming a maximum 
fractional energy of the virtual photon of y ma x = 0.95. 
We only consider c.m.s. energies which allow one to ac- 
cess x values at least down to 10 -3 even for a minimum 
Q 2 = 2.5 GeV 2 to achieve the goal of constraining helicity 
PDFs in the so far unexplored small x region. 



TABLE I: Combinations of electron and proton energies used 
in our analyses and the corresponding c.m.s. energies and min- 
imum values of x accessible for Q 2 — 1 and 2.5 GeV 2 and 
2/max = 0.95. For each data set a modest integrated luminos- 
ity of 10fb _1 is assumed. 
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FIG. 1: [color online] Kinematic range in x and Q 2 accessible 
with two different c.m.s. energies at an EIC for 0.01 < y < 
0.95 (hatched areas). The points illustrate the coverage of 
currently available data from longitudinally polarized DIS and 
pp experiments (see text). 

Figure [T]illustrates the kinematic coverage of an EIC in 
the x — Q 2 plane for both a first stage with yfs ~ 45 GeV 
and a full energy eRHIC option with y/s ~ 140 GeV. An 
EIC will dramatically extend both the x and the Q 2 cov- 
erage of existing fixed-target DIS experiments, denoted 
by the different symbols, by about two orders of magni- 
tude, opening up unprecedented opportunities for preci- 
sion studies of helicity PDFs. Present QCD analyses [l|- 
IH need to include all DIS data down to Q 2 = 1 GeV in 
order to reach x values of about 4.6 x 10 -3 . Allowing for a 
higher, more conservative cut in Q 2 , say Q 2 > 2.5 GeV 2 , 
to study the possible relevance of "higher twist" cor- 
rections, which are suppressed by inverse powers of the 
hard scale Q, search for possible deviations from standard 
DGLAP scale evolution [13j, or to test the applicability of 
the assumed factorized pQCD framework, is not possible 



for the time being as it would limit the accessible x range 
too much. As can be inferred from Tab. [J and Fig. [TJ at 
an EIC one can perform all these important studies and 
can still reach down to at least x ~ 0.001 even for the 
lowest c.m.s. energy option we are considering here. 

Also shown in Fig.Q]are the data from polarized pp col- 
lisions at BNL-RHIC [9J which currently provide the best 
constraint on the polarized gluon distribution. We note 
that assigning a transverse momentum dependent observ- 
able in pp collisions to a single value of x probed in the 
gluon is a gross oversimplification. Choosing x — 2px / \fs 
nevertheless gives at least a rough idea about the lowest 
x values accessible in pp collisions for a given px ■ Also in 
DIS at NLO accuracy, information on PDFs is contained 
within convolutions with hard scattering cross sections, 
and x merely reflects the lowest possible momentum frac- 
tion accessible in experiment. The lever-arm in Q 2 for 
any given value of x at an EIC will allow one to ex- 
tract information on Ag(x,Q 2 ) from scaling violations, 
i.e., the rate at which the polarized DIS cross section 
changes with scale Q for a given fixed value of x. De- 
pending on the lowest value of Q 2 used in the analyses, 
such studies can be pushed down to about i~lx 10~ 4 
as will be demonstrated below. 

Another unique opportunity at an EIC, although not 
being pursued further in this study, is to access novel 
electroweak structure functions in DIS which probe com- 
binations of quark helicity PDFs at medium-to-large x 
and large Q 2 ~ M{y different from those constrained by 
DIS data at much lower scales and described solely by 
the exchange of a virtual photon. Clearly, the kinematic 
coverage along with the envisioned unprecedented lumi- 
nosity and the possibility of having polarized beams are 
the biggest assets of an EIC which make a suite of preci- 
sion QCD studies possible [7|, ranging from the proton's 
helicity structure considered in this paper to hadron "to- 
mography" through exclusive processes. 

To quantify the impact of an EIC on our understanding 
of helicity PDFs, we generate sets of pseudo-data for each 
of the three c.m.s. energies listed in Tab. [TJ We use the 
PEPSI Monte Carlo (MC) generator [27[ to produce ficti- 
tious EIC data for the inclusive and semi-inclusive DIS of 
longitudinally polarized electrons and protons with iden- 
tified charged pions and kaons in the final-state. We 
demand a minimum Q 2 of 1 GeV 2 , a squared invariant 
mass of the virtual photon-proton system larger than 
W 2 = 10 GeV 2 , and 0.01 < y < 0.95. The range of y 
is further restricted from below by constraining the de- 
polarization factor of the virtual photon 



D{y) 



y(y-2) 



2(l-y)(l + R) 



y 



(2) 



to be larger than 0.1. R denotes the ratio of the lon- 
gitudinal to transverse virtual photon cross sections. To 
ensure detection of the scattered lepton we require a min- 
imum momentum of 0.5 GeV, and, in case of SIDIS, only 
hadrons with a momentum larger than 1 GeV and a frac- 
tional energy in the range 0.2 < z < 0.9 are accepted. 
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All particles detected in the final-state should be at least 
1 degree away from the beam directions. The statistical 
accuracy of each DIS and SIDIS data set corresponds to 
a modest accumulated integrated luminosity of 10fb _1 , 
equivalent to about one to two months of operations for 
the anticipated luminosities for cRHIC 26(, except for 
the 5 x 100 GeV option which requires about a year of 
running. 

The PEPSI MC '27} is currently the only code that 
allows one to generate events with definite helicities of 
the colliding lepton and proton beams, i.e., to study the 
longitudinal double-spin asymmetry 



A l{ (x,Q 2 ) = 



da++ - da+~ 
da++ + da+- 

D(y) 9l{x > Q2) 



F 1 {x,Q 2 ) 



(3) 
(4) 



which is related to the ratio of virtual photoabsorption 
cross sections, expressed by DIS structure functions in 
(|4]), through the depolarization factor D{y). In (j4j, and 
also in @, we have neglected kinematic corrections pro- 
portional to 7 = y / 4M 2 x 2 /Q 2 , with M the proton mass, 
which are negligible at a collider as can be deduced from 
the x — Q 2 coverage depicted Fig. [1] While containing 
spin-dependent hard scattering matrix elements at LO 
accuracy, the PEPSI MC is not capable of simulating 
parton showers which properly track the polarizations 
of the partons involved, and hence this option has been 
turned off for generating the EIC data. QED radiative 
corrections are known to be sizable and complicate the 
determination of the "true" values of x and Q 2 . On the 
other hand, we have learned a great deal on how to con- 
trol and unfold these corrections from years of successful 
DIS structure function measurements at DESY-HERA. 
Undoubtedly, available MC tools [28j will be further re- 
fined in the upcoming years, and we do not consider QED 
radiative corrections to be a major limitation on pro- 
posed DIS and SIDIS measurements at an EIC. We note 
that the typical size of the double spin asymmetry (J3]) at 
the lowest x values accessible at an EIC can be as small 
as a few times 10 -4 , depending on the yet unknown be- 
havior of Ag(x, Q 2 ) in this kinematic regime. This size 
sets the scale at which one needs to control systematic 
uncertainties due to detector performance or luminos- 
ity measurements. Most likely, the dominant source of 
systematic uncertainty will be the determination of the 
beam polarizations which will lead to a scale uncertainty 
in spin asymmetry measurements. In any case, DIS mea- 
surements, even with the anticipated high precision, will 
be far from being the most challenging measurements to 
be performed at an EIC Q. 

Monte Carlo data for the ratio gx/F\ in DIS and SIDIS 
are generated in 4 [5] bins per decade in Q 2 [x] spaced 
logarithmically. As the actual pseudo-data used in our 
global analyses, we take the ratio Pi/Fi computed at 
NLO accuracy using the DSSV+ [l| and MRST [U 
polarized and unpolarized PDFs, respectively, and as- 



sign to each (x, <5 2 )-bin the same relative statistical un- 
certainties as obtained with the MC event generator cor- 
responding to an integrated luminosity of 10fb _1 and 
assuming 70% beam polarizations. In addition, we ran- 
domize the pseudo-data in each bin within these one- 
sigma uncertainties. In total we add 234 data points 
for DIS and about 800 points for SIDIS to the existing 
DSSV/DSSV+ global analysis framework 0, Q]| based 
on 570 DIS, SIDIS, and pp data. For the SIDIS data 
with identified charged pions and kaons we assign an ad- 
ditional, conservative 5 and 10% relative uncertainty to 
the EIC pseudo-data to reflect our current incomplete 
knowledge of parton-to-pion and parton-to-kaon frag- 
mentation functions, respectively, based on uncertainty 
estimates for the DSS sets of FFs in Ref. [3(| • As already 
mentioned, various upcoming data sets are expected to 
greatly advance our knowledge of FFs in upcoming years 
though. If necessary, this can be supplemented with 
measurements of unpolarized hadron multiplicities at the 
EIC. Some expectations for charged kaon production in 
unpolarized SIDIS can be found in Q. 



To quantify the impact of the simulated EIC data on 
our understanding of the spin structure of the nucleon, we 
first need to define some set of reference results which re- 
flects our current knowledge of helicity PDFs, including a 
faithful estimate of their present uncertainties. Here, we 
follow the framework of the original DSSV global anal- 
ysis [l| and keep the same functional form, initial scale 
Qo = 1 GeV, unpolarized reference PDFs, standard \ 2 
minimization procedure, and number of free fit param- 
eters to facilitate comparisons. However, as mentioned 
above, we update the data sets by including COMPASS 
DIS [3l| and SIDIS [l6[ data which became available only 
after the DSSV analysis was completed. The resulting 
new best fit, labeled as DSSV+, which hardly differs from 
the published DSSV results, is used as our baseline fit. 
Despite ongoing experimental efforts, the obtained helic- 
ity PDFs will reflect to a good approximation of what will 
be known by the time an EIC will start operating. Some 
improvements are expected in the meantime at large x 
and low Q 2 from JLab-12 and towards somewhat smaller 
values of x than indicated in Fig. [1] from ongoing polar- 
ized RHIC pp running at y/s — 500 GeV. Once a suf- 
ficient amount of data has been accumulated, measure- 
ments of single-spin asymmetries in W boson production 
at RHIC are expected to improve u and d quark and an- 
tiquark helicity PDFs at Q ~ Mw and medium-to-large 
values of x [32] . As in Ref. [l[ , we will use primarily the 
robust Lagrange multiplier method [24| to quantify un- 
certainties with and without including the simulated DIS 
and SIDIS data. In addition, once EIC data are included 
in the fit, the standard Hessian method (25j, which ex- 
plores the vicinity of the x 2 minimum in the quadratic 
approximation, also starts to produce reliable results and 
can be compared to the results obtained with Lagrange 
multipliers. 
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FIG. 2: [color online] Projected EIC data for the structure 
function gi(x, Q 2 ) for three different combinations of electron 
and proton energies. Constants are added to gi to separate 
the different x bins. The solid lines are the result of the 
DSSV+ best fit, and the shaded bands illustrate the current 
uncertainty estimate. Multiple data points at a given x, Q 2 
are displaced horizontally to make them more easily visible. 
The hatched triangular area indicates the region covered by 
present data. 



III. IMPACT OF DIS AND SIDIS DATA 

Figure[5]illustrates our simulated data sets for inclusive 
polarized DIS at an EIC for the three different choices 
of c.m.s. energies listed in Tab. HI The error bars were 
determined as outlined in the previous Section and re- 
flect the expected statistical accuracy for a modest inte- 
grated luminosity of 10 fb _1 . As indicated by the hatched 
area, existing fixed target DIS data (see Fig. [1]) populate 
only the lower left corner of the kinematic plane but con- 
nect well or overlap with the lowest Q 2 values accessible 
with the 5 x 100 GeV data set. Relaxing our conservative 
constraint on the depolarization factor ([2]), D(y) > 0.1, 
would significantly increase the overlap to even lower val- 
ues of Q 2 . We note in passing that if one can control 
systematic uncertainties very well at an EIC, which is 
definitely the goal, one might try to aim for polarized 
cross section rather than asymmetry measurements in 
the future. This would have the added benefit of being 



independent of the ratio R of the longitudinal to trans- 
verse virtual photon cross sections. The shaded bands 
in Fig. [2] correspond to the current uncertainties as esti- 
mated in the DSSV analysis based on the Lagrange multi- 
plier method. At low x, outside the range constrained by 
present data, these bands essentially reflect the flexibility 
of the chosen functional form and are a mere extrapola- 
tion. 

As is already obvious from Tab. HI DIS measurements 
for 20 x 250 GeV collisions are crucial to reach x values of 
around 10 -4 while still maintaining at least some lever- 
arm in Q 2 . With energies of up to 5 x 250 GeV, envisioned 
in the first stage of eRHIC, one can still cover x values 
down to 5 x 10~ 4 for Q 2 > 2.5 GeV 2 . Having available 
an as large as possible range in Q 2 for any given fixed 
value of x is of outmost importance for studying scal- 
ing violations which are a key prediction of pQCD. Even 
though the DIS structure function 171 probes mainly the 
sum of quark and antiquark PDFs, its scaling violations 
at small enough values of x are approximately related to 
the polarized gluon density, 



dgi(x, Q 2 
d\nQ 2 



-Ag(x,Q 2 



(5) 



which underlines the importance of precisely measuring 
them. In very much the same way, unpolarized DIS data 
from the DESY-HERA experiments HI and ZEUS [H 
provide the best constraint on the gluon density at small 
momentum fractions in all global QCD analyses thanks 
to their vast range in x and Q 2 only accessible at collider 
energies. It is fair to say, that with presently available 
polarized DIS data one can hardly utilize the relation 
(O to determine A<? because of the much too limited 
kinematic coverage. 

The prospects for measuring dgi(x,Q 2 )/dlnQ 2 at an 
EIC are summarized in Fig. [3] The projected scaling vi- 
olations are obtained from the DIS pseudo-data shown 
in Fig. [2] For a given bin in x, one needs, of course, at 
least measurements of gi(x,Q 2 ) at two different values 
of Q 2 which are precise enough to reliably determine the 
derivative dgi(x, Q 2 )/d\nQ 2 from a difference quotient. 
For the binning in x and Q 2 adopted in our analysis and 
the assumed integrated luminosity of 10fb _1 , a measure- 
ment of dgi(x, Q 2 )/dhiQ 2 down to about x ~ 1 x 10~ 4 
appears to be conceivable assuming 20 x 250 GeV col- 
lisions. Likewise, a first stage option of an EIC with 
5 x 250 GeV will have sensitivity down to x ~ 5 x 10~ 4 . 
This also roughly delineates the range in x where one 
can expect to put a sensible constraint on Ag(x, Q 2 ) with 
an inclusive DIS measurement at an EIC assuming that 
([5]) is a good approximation. The smallness of the pro- 
jected statistical errors indicates that all inclusive and 
semi-inclusive DIS measurements discussed here are sys- 
tematics limited. Precision measurements will require 
a percent-level control of the many different sources of 
systematic uncertainties such as the luminosity and po- 
larization measurements but also of the resolution and 
calibration of the required detector elements and in the 
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FIG. 3: [color online] Theoretical expectations for the logarithmic scaling violations dgi(x,Q 2 )/dlogQ 2 based on the DSSV+ 
best fit (solid lines) in different bins of x. The shaded bands reflect the current (asymmetric) uncertainties of the DSSV helicity 
PDFs. The points illustrate the expected accuracy for measurements at an EIC with 5 GeV electrons (triangles) based on the 
projected data shown in Fig. [2] The open circles include also results obtained with 20 GeV electrons. Data points at the same 
given (x,Q 2 ) are displaced horizontally to distinguish them better. 



unfolding of QED radiative corrections. 

Our projected SIDIS data for identified charged pions 
and kaons share the same x and Q 2 binning as the DIS 
data presented in Fig. [3] but have slightly larger uncer- 
tainties since we assign an up to 10% additional rela- 
tive uncertainty due to FFs as explained in Sec. II. The 
SIDIS cross section for a hadron h can be expressed by a 
structure function g±(x, z,Q 2 ), h — 7r ± , -ftT ± , which now 
depends also on the fraction z of the momentum of the 
fragmenting quark or gluon taken by the observed hadron 
h. gi(x,z,Q 2 ) exhibits similar scaling violations as the 
inclusive g\ and will contribute to constraining Ag in a 
global QCD analysis but SIDIS data draw their relevance 
from their sensitivity to different quark and antiquark 
flavors. To calibrate measurements of polarized SIDIS at 
an EIC, one will look at unpolarized hadron multiplicities 
first, which on the one hand will help to map out down to 
which values of Q 2 the leading-twist pQCD framework is 
a good approximation for SIDIS and on the other hand 
will test and improve our knowledge of FFs. 

Figure @] demonstrates the impact of the projected 
combined EIC data for DIS and SIDIS on extractions 
of the polarized sea quark and gluon densities. Here we 
utilize only data which can be obtained already with the 
initial stage of the eRHIC option, i.e., for collision ener- 
gies of 5 x 100 and 5 x 250 GeV. The outer bands in each 
panel refer to the present ambiguities for helicity PDFs 
as determined in the DSSV analysis 0, d, [l5| and cor- 
responds to a conservative increase in the total x 2 used 
to determine the goodness of the fit by nine units. This 
value of Ax 2 was regarded to lead to a faithful estimate 
of present uncertainties in (ll 13, [lH . The smaller, inner 
bands are obtained with the same global analysis frame- 
work, functional form for the PDFs, number of free fit 
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FIG. 4: [color online] Impact of projected EIC data for the 
DIS and SIDIS of 5 GeV electrons on 100 and 250 GeV pro- 
tons on the determination of helicity sea quark PDFs and the 
gluon. The light shaded bands illustrate present uncertainty 
estimates and the dark shaded, inner bands the improvements 
expected due to EIC data (see text). 

parameters, and A% 2 criterion but now include also the 
projected EIC data. As can be seen, the expected im- 
provements are dramatic, in particular, for the polarized 
gluon density below x ~ 0.01 but also for the individual 
sea quark flavors. We wish to emphasize that as long 
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FIG. 5: [color online] x 2 profiles for the first moments of he- 
licity sea quark PDFs and the gluon truncated to the region 
0.001 < x < 1. The results are based on using only cur- 
rent data (DSSV+) and sets of projected EIC data with two 
different c.m.s. energies. 



as we limit ourselves to the range x > 10 -3 , the results 
shown in Fig. |4] do not require to analyze data below 
Q 2 ~ 2.5 GeV where the perturbative framework eventu- 
ally starts to become unreliable and/or where 1/Q sup- 
pressed power corrections may become relevant. At an 
EIC one can systematically study the validity of the lead- 
ing twist pQCD framework assumed in all global QCD 
analyses so far by varying the lower cut-off scale Q m i n 
above which one starts to include data in the fit. It 
should also be stressed that only the relative improve- 
ment of the uncertainties in Fig. 0] i.e., the differences 
between the inner and outer error bands, is of signifi- 
cance here for estimating the physics impact of an EIC 
since the generation of the pseudo-data requires to as- 
sume a certain set of polarized PDFs. Of course, only 
real EIC data will eventually reveal the actual functional 
form of the helicity PDFs at small x. 

Figure [5] visualizes the improvements due to EIC data 
in terms of truncated moments of helicity PDFs, which 
are also used to determine the uncertainty bands in Fig.j4] 
with the help of the Lagrange multiplier method 
ITsl l2~ij . For each parton flavor / we minimize x 2 with an 
additional constraint on the value of its truncated first 
moment 

Af{Q 2 ,x min ,x ma ^) = / Af(x,Q 2 )dx (6) 

Smin 

implemented through a Lagrange multiplier. In this 
way we can map out the x 2 profile as a function of 
Af{Q 2 , x m in, Xmax) away from its best fit value without 




0.018 0.02 0.022 0.024 -0.04 -0.03 -0.02 -0.01 
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FIG. 6: [color online] As in Fig. [5j but now evaluated in the 
range 0.0001 < x < 0.01 and using also a set of projected EIC 
data for collisions of 20 GeV electrons on 250 GeV protons. 



any restrictions on the parameter space. Thereby we gen- 
erate a large set of alternative PDFs for each point along 
the x 2 contour. Upon choosing a certain maximum in- 
crease A% 2 , which is still tolerated for a good fit, one 
arrives at the uncertainty bands shown in Fig. 2) 

A more direct way of estimating PDF uncertainties is 
the standard Hessian method [25j which is based on the 
assumption that the \ 2 profiles are quadratic in the vicin- 
ity of their minima. As can be inferred from the profiles 
in Fig. [5j the truncated moments A/(Q 2 , 0.001, 1) are 
only very weakly constrained by presently available data, 
and uncertainties are very large. Clearly, the quadratic 
approximation does only work well for not too large A\ 2 
[2j, and, hence, reliable Hessian eigenvector PDF sets 
for Ax 2 = 9 cannot be constructed to estimate uncer- 
tainties. However, including just one of the projected 
EIC data sets not only considerably reduces uncertain- 
ties for A/(Q 2 , 0.001, 1), which can be conveniently read 
off from the width of the x 2 profiles at any desired value 
of A% 2 , but also leads to approximately parabolic x 2 
profiles. Hence, to check the consistency of our error 
estimations, we also determined the PDF uncertainties 
with the now applicable Hessian method by construct- 
ing appropriate eigenvector PDF sets corresponding to 
Ax 2 = 9. We find very similar, basically undistinguish- 
able results as for the inner uncertainty bands shown in 
Fig. |4] and obtained with Lagrange multipliers. 

It also turns out, see Fig. [SJ that helicity PDFs are 
already well constrained down to x = 1 X 10~ 3 by EIC 
data for 5 x 100 GeV collision because essentially all x- 
bins at that particular c.m.s. energy fall into the region 
x 10~ 3 . Nevertheless, additional data for 5 x 250 GeV 
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collisions will further improve the constraint on Ag, 
mainly because of the extended lever-arm in Q 2 for stud- 
ies of scaling violations, see Fig. [3j 

To explore the impact of projected EIC data for colli- 
sions of 20 GeV electrons on 250 GeV protons, envisioned 
at a full energy eRHIC, we perform a similar analysis as 
in Fig. [5] but now for the moments truncated in the range 
from x = 0.0001 to 0.01. As has to be expected, current 
constraints are even weaker than for the range x > 0.001 
considered in Fig. [5j resulting in essentially flat \ 2 P ro ~ 
files. Again, projected EIC data will lead to dramatic 
improvements even at an initial stage with only 5 GeV 
electrons available. In particular for precision studies of 
Ag at low x, the 20 x 250 GeV data will be extremely 
crucial as they greatly enlarge the x range where scaling 
violations can be studied, see Fig. [3J This is reflected by 
the significant further reduction of the uncertainties in 
the lower right panel of Fig. [5] In case that the helic- 
ity PDFs exhibit some sign change at medium-to-large 
x, as, for instance, Au(x Q 2 \ As(x, Q 2 ), and Ag(x,Q 2 ) 
in the DSSV analysis [l|, @, LL5J, a numerically signifi- 
cant contribution to their first moments may arise from 
the small x region, i.e., x < 0.01, only accessible at an 
EIC. Having 20 x 250 GeV data at hand, one can extend 
the uncertainty bands shown in Fig. [3] to x = 1 x 10~ 4 
and perhaps to even lower values of x if data down to 
Q 2 ~ 1 4- 1.5 GeV 2 appear to be amenable to standard 
leading twist factorization and pQCD methods. 

The much reduced uncertainties of helicity PDFs 
thanks to EIC data will allow one to quantitatively ad- 
dress most of the physics questions concerning the pro- 
ton's spin structure raised in the Introduction. From the 
X profiles shown above, one can gather that one can 
determine the small x behavior of As (a;, Q 2 ) very accu- 
rately, mainly from SIDIS kaon data. In the DSSV anal- 
ysis U 0, EH , As acquires most of its sizable negative 
x-integral in the so far unmeasured small x region in or- 
der to respect a constraint from SU(3) flavor symmetry. 
The latter is expressed in terms of the two hyperon decay 
matrix elements F and D which are experimentally well 
known, explaining the rather small uncertainty band for 
As(x, Q 2 ) at small x. However, the applicability of this 
constraint has been questioned, and large SU(3) break- 
ing effects are certainly not excluded yet (l7l]la|. The 
ambiguities related to the assumptions about the small x 
behavior of As, i.e., the amount of SU(3) breaking, also 
drives the current uncertainties of the first moment of the 
flavor singlet combination AS [l], 0, EH which enters in 
the proton spin sum rule. 

From the \ 2 profile in the lower left panel of 
Fig. [6] one can read off that the truncated moment 
As(Q 2 , 0.0001, 0.01) is indeed sizable and negative in the 
DSSV analysis. EIC data are expected to constrain it to 
within about 5%. Clearly, even with modest systematic 
uncertainties present, one can easily quantify the amount 
of SU(3) flavor breaking effects at an EIC which is a 
crucial ingredient to our understanding of the partonic 
structure of hadrons and the possible relevance of chiral 



corrections as estimated, for instance, in the cloudy bag 
model [H. 

Another interesting question related to strangeness is 
a possible asymmetry As(x, Q 2 ) — As(x, Q 2 ) ^ which 
is also one of the least well determined quantities in case 
of unpolarized PDFs. At an EIC, the difference between 
yields for K + and K~~ will provide sensitivity to such 
kind of quantities but likely requires an improved un- 
derstanding of the analyzing power given by the ratio 
of the favored to unfavored strangeness fragmentation 
functions and , respectively. A first feasibility 

study for unpolarized SIDIS can be found in Q but due 
to the current limitations for FFs we do not pursue this 
further here. The LHC is already starting to provide in- 
teresting new insights into unpolarized strangeness dis- 
tributions from precision measurements of electroweak 
boson production [34| which can be utilized at the EIC 
to first improve our knowledge of kaon FFs in unpolar- 
ized SIDIS. This information should be then sufficient to 
study As{x, Q 2 ) - As(x, Q 2 ). 
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FIG. 7: [color online] As in Fig.[4]but now for the difference of 
the light sea quark densities Au — Ad. The dotted line shows 
an expectation from the chiral quark soliton model and 
the dashed line the corresponding asymmetry for unpolarized 
PDFs from the CTEQ6 analysis [13] multiplied by (-1). 

The asymmetry in the light quark sea, Au(x,Q 2 ) — 
Ad(x,Q 2 ), is of particular interest as well. Firstly, it 
is know to be sizable in the unpolarized case [23| and 
secondly it can be predicted in various models of the nu- 
cleon structure such as the chiral quark soliton model 
[22| where one expects an SU(2) flavor breaking of the 
sea which is at least as large as what has been already 
observed for unpolarized PDFs but with its sign reversed. 
Figure [7] shows both an expectation from the chiral quark 
soliton model [22J and a ty pica l breaking obtained in un- 
polarized global PDF fits [23| . As in Fig. |U the larger 
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(outer) error band corresponds to an uncertainty esti- 
mate for Au(x,Q 2 ) — Ad(x,Q 2 ) by DSSV based on a 
Lagrange multiplier analysis of currently available data 
and a tolerated increase in \ 2 by 9 units. The impact 
of projected EIC data for 5 x 100 and 5 x 250 GeV colli- 
sions is illustrated by the inner error band. This exercise 
shows that asymmetries Au(x, Q 2 ) — Ad(x, Q 2 ) of about 
0.02 can be resolved, which is more than sufficient to test 
typical model expectations. We recall that even in the 
absence of a non-perturbative asymmetry at some low 
hadronic scale, i.e., Am = Ad, a non-zero asymmetry 
will be generated perturbatively through QCD scale evo- 
lution at NLO accuracy [35| and beyond. Likewise, at 
NNLO an x dependent, local s(x, Q 2 ) — s(x,Q 2 ) asym- 
metry will develop under QCD scale evolution [361 ] . 
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FIG. 8: [color online] Correlated uncertainties for the first 
moments of the flavor singlet combination AE and the gluon 
helicity density Ag computed in the region 0.001 < x < 1. 
The green, red, and yellow shaded areas are based on fits 
to current data and to projected EIC data with 5 GeV and 
20 GeV electron beams, respectively. The symbol denotes the 
DSSV+ best fit. 

Finally, we look into what can be achieved for the 
first moments of the flavor singlet combination AS and 
the gluon helicity density Ag which both enter the pro- 
ton spin rule. Figure [8] shows the correlated uncertain- 
ties for the truncated moments computed in the region 
0.001 < x < 1 with and without including projected EIC 
data sets. As for all our studies presented above, an EIC 
will greatly improve the uncertainties, in particular, for 
Ag which is essentially unconstrained so far. As can be 
seen, Ag(Q 2 , 0.001, 1) and AS(Q 2 , 0.001, 1) can be con- 
strained up to about ±0.05 and ±0.02, respectively, if 
20 x 250 GeV data are included in the PDF analyses. 
However, already at the initial stage of an EIC a very 
significant reduction of uncertainties can be achieved. 



Again, as a check, very similar uncertainty estimates 
have been obtained with the Hessian method once pro- 
jected EIC data are included in the global analysis frame- 
work. Similar improvements as in Fig. [5] are found 
for the truncated moments for x > 1 x 10 -4 but then 
20 x 250 GeV data are essential, cf. the results for Ag 
shown in Fig. [51 Although already mentioned in the In- 
troduction, we recall that the decomposition of the pro- 
ton spin in Quantum Field Theory is non-trivial and still 
under active discussion While the flavor singlet con- 
tribution AE appears universally in all proposed spin 
sum rules, the first moment of the helicity gluon density 
acquires the interpretation as the gluon spin contribution 
to the proton spin only in light-cone gauge which is the 
natural gauge to define parton distributions. Quark and 
gluon orbital angular momenta will be part of another 
suite of unique measurements at an EIC aiming at the 
nucleons spatial structure @- 



IV. OTHER OPPORTUNITIES IN DIS 

The presented examples for a physics program with 
longitudinally polarized lepton-nucleon collisions at an 
EIC are just the tip of the iceberg and only cover the most 
important measurements which can be done in DIS and 
SIDIS. There are many other avenues at an EIC which 
can be pursued to further our understanding of the he- 
licity structure of the nucleon. For instance, at high Q 2 
an EIC gives unique access to electroweak effects in po- 
larized DIS which have not been measured so far. Pre- 
liminary studies can be found in [7] . Such measurements 
can be also extended to SIDIS. In general, electroweak 
DIS results will constrain helicity PDFs at medium-to- 
large values of x but at much higher scales Q 2 ~ 
than currently existing fixed-target DIS data which can 
be described solely by one-photon exchange. Although 
the QCD scale evolution is expected to work well at large 
Q 2 , it has not been tested yet for helicity PDFs. More 
importantly, such measurements should provide a clean 
way to quantify possible higher twist contributions at 
large momentum fractions x and low Q 2 from combined 
fits. 

Another important observable is charm production in 
polarized DIS, i.e., the charm contribution g1(x,Q 2 ) to 
the inclusive structure function g\. So far a proper- 
treatment of heavy flavors in polarized DIS is irrelevant 
since in the presently covered (x, Q 2 ) range its contri- 
bution to gi amounts to less than 1%. At smaller x 
values, accessible for the first time at an EIC, the size 
of gf(x, Q 2 )/ gi(x, Q 2 ) very much depends on what is as- 
sumed for the helicity gluon density as charm is domi- 
nantly produced through photon-gluon fusion j*g — » cc. 
For a small Ag, as, for instance, in the best fit of DSSV, 
g1(x,Q 2 ) remains to be negligible but can account for 
as much as 15% of the inclusive g\ at x = 0.001 for a 
larger gluon distribution; see Q for some quantitative 
estimates. More theoretical work is clearly needed here, 
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since the relevant cross section for 7*5 — > cc with full de- 
pendence on the charm quark mass has been calculated so 
far only at LO accuracy [37j , Also, variable flavor num- 
ber schemes have not been considered yet for polarized 
PDF sets and need to be developed in the future. 

It is also conceivable that an EIC can store polarized 
Helium-3 beams which essentially act as source of polar- 
ized neutrons as | 3 He) = 0.865|n) + 2(-0.027)|p) and if 
the spectator protons in an inelastic collision are detected 
with the help of Roman pots. The prime physics motiva- 
tion for studying longitudinally polarized lepton-neutron 
collisions is not so much the extra handle on the flavor 
decomposition but mainly the Bjorken sum rule [38] 

J dx [<^(x,Q 2 - 5l "(x,g 2 )] = l -C Bj [a s {Q 2 )]g A (7) 

which is currently experimentally verified to about 10%. 
The Bjorken sum rule is a rare example of a well- 
understood quantity in pQCD with O(af) corrections 
to Cb 3 [a s (Q 2 )] being known [39| and potentially large 
1/Q 2 higher twist corrections expected to be small in 
the perturbative regime [4(|. Of course, significantly im- 
proving the current level of experimental accuracy re- 
quires, among other things, percent level control for 3 He 
polarimetry which needs some novel technical ideas. In 
Ref. it has been estimated that one has to access x val- 
ues down to 0(1O~ 4 ) to limit extrapolation uncertainties 
in the non-singlet combination Aq 3 = Au + Au — (Ad + 
Ad) effectively probed by the Bjorken sum rule (0 to 
a level of 1 — 2%. Further theoretical interest in the 
Bjorken sum is generated from its relation to the Adler 
function in e + e~ annihilation through the Crewther rela- 
tion [4l[ which has been worked out up to 0(af) recently 
[3^ | . In terms of providing novel information on the fla- 
vor separation of helicity PDFs, neutron data may help 
to reduce uncertainties for Ad and Ad beyond what can 
be achieved with SIDIS data in polarized electron-proton 
collisions thanks to the u f-> d isospin rotation. If tech- 
nical issues concerning 3 He polarimetry do not prove to 
be too demanding, a quantitative estimate of the impact 
3 He data on the determination of helicity PDFs can be 
made along very similar lines as in our study. Statistical 
uncertainties will be comparable to the once obtained for 
polarized protons in Fig. [2J 

Finally, precision QCD studies of the helicity struc- 
ture of nucleons at an EIC may reveal tensions with 
DGLAP scale evolution [l3j which are expected at suffi- 
ciently small but otherwise hard to pinpoint values of x 
42]. In contrast to the unpolarized case, the dominant 
contribution of gluons mixes with quarks also at x <C 1. 
From the standard scale evolution |l3[ one expects for 
the small x behavior 

Aq(x,Q 2 )Ag(x,Q 2 )^ 

exp [const x a s ln(Q 2 //i 2 ) ln(l/x)] 1/2 (8) 

assuming for simplicity a fixed coupling a s . In [42} it 
was demonstrated that this simple behavior can strongly 



underestimate the rise at small x due to other poten- 
tially large double logarithmic contributions of the type 
a s ln 2 (l/x) n in the n-th order of a s which are beyond the 
standard framework. This gives rise to a power-like be- 
havior of c/i at small x of the form gi (x, Q 2 ) ~ (l/x)°^ as \ 
There are qualitative arguments that in the polarized 
case these logarithms in 1/x are more relevant than in 
the unpolarized case [42j ■ Only data can eventually re- 
veal if the kinematic reach of an EIC is sufficient to actu- 
ally observe deviations from conventional scale evolution 
in polarized DIS. 



V. SUMMARY 

We have presented a detailed quantitative assessment 
of the dramatic impact a future EIC will have on deter- 
minations of helicity quark and gluon densities and their 
contributions to the proton spin. Key asset of a first po- 
larized lepton-nucleon collider will be its unprecedented 
kinematic coverage both down to small momentum frac- 
tions x~lxl0~ 4 and to large scales Q, implying a suffi- 
ciently large c.m.s. energy of the collisions. This is essen- 
tial to further our understanding of the nucleon's helicity 
structure to level which is sufficient to quantitatively ad- 
dress outstanding questions about the role of polarized 
gluons and the flavor structure of sea quark densities at 
small momentum fractions x. The necessary precision 
measurements in polarized inclusive and semi-inclusive 
deep-inelastic scattering only require modest integrated 
luminosities but good control over all sources of system- 
atic uncertainties ranging from luminosity and polariza- 
tion measurements, detector acceptance and resolution, 
to a proper unfolding of QED radiative corrections. 

All presented results were obtained by performing a 
series of global QCD analyses at NLO accuracy based on 
realistic sets of pseudo-data for the inclusive and semi- 
inclusive deep-inelastic scattering of longitudinally po- 
larized electrons and protons at different, conceivable 
center-of-mass system energies. The dramatic physics 
impact of such data sets has been quantified by esti- 
mating uncertainties for all relevant quantities with the 
robust Lagrange multiplier method and by comparing 
them to present-day helicity PDF uncertainties. An EIC 
will provide precise information on the helicity depen- 
dent gluon and flavor separated quark densities down the 
momentum fractions of about 10 -4 which in turn will ac- 
curately determine their contribution to the spin of the 
proton. 

We have briefly highlighted other interesting opportu- 
nities related to helicity PDFs which can be only pursued 
at an EIC such as a precision measurement of the Bjorken 
sum role which requires, however, to overcome all techni- 
cal challenges related to the need for having an effective 
polarized neutron beam. Charm and electroweak con- 
tributions to helicity DIS structure functions are other 
prominent examples for measurements uniquely tied to a 
high energy polarized lepton-nucleon collider. 
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